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Abstract 

We present a study of the neutron electric dipole moment (djv) within the framework of lattice 
QCD with two flavors of dynamical light quarks. The dipole moment is sensitive to the topological 
structure of the gauge fields, and accuracy can only be achieved by using dynamical, or sea quark, 
calculations. However, the topological charge evolves slowly in these calculations, leading to a rel¬ 
atively large uncertainty in djsf. It is shown, using quenched configurations, that a better sampling 
of the charge distribution reduces this problem, but because the CP even part of the fermion de¬ 
terminant is absent, both the topological charge distribution and djsf are pathological in the chiral 
limit. We discuss the statistical and systematic uncertainties arising from the topological charge 
distribution and unphysical size of the quark mass in our calculations and prospects for eliminating 
them. 

Our calculations employ the RBC collaboration two flavor domain wall fermion and DBW2 gauge 
action lattices with inverse lattice spacing a~^ ~ 1.7 GeV, physical volume F ~ (2 fm)^, and light 
quark mass roughly equal to the strange quark mass (m^ea = 0.03 and 0.04). We determine a value 
of the electric dipole moment that is zero within (statistical) errors, |dAr| = —0.04(20) e-O-im at 
the smaller sea quark mass. Satisfactory results for the magnetic and electric form factors of the 
proton and neutron are also obtained and presented. 

PACS numbers: 11.15.Ha, ll.30.Er, ll.30.Rd, 12.38.Aw, 12.38.-t 12.38.Gc 
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I. INTRODUCTION 


One of the most intriguing aspects of quantum chromodynamics (QCD) is that it allows 
a gauge invariant interaction term that is separately odd under time-reversal (T) and parity 
(P) transformations, the so-called 6 term. The presence of such a term has the profound 
effect that the Strong interactions violate the combined symmetry charge-conjugation (C) 
times P. The existence of P and T violating interactions in the action imply permanent 
electric dipole moments for fundamental particles. Presently, the most precise search for 
a permanent electric dipole moment comes from the measurement of the electric dipole 
moment of the neutron, dj^j. In the Standard Model, the CP-odd phase of the CKM mixing 
matrix also produces a non-vanishing value for d^, but only beyond one loop order in 
the Weak interaction. Consequently, this contribution to is estimated to be less than 


10“^° e-cm, many orders of magnitude below the current experimental bound|l|, djsi = \disf \ < 
6.3 X 10“^® e-cm (see also j^). There are recent proposals to improve this bound by two 
to three orders of magnitude by studying the electric dipole of the deuteron at Brookhaven 
National Laboratoryjsf and an isotope of radium (^^®Ra) at Argonne National Laboratory 
The latter is now underway. 

Using this experimental bound with theoretical estimates of dN/e BSflBSQ, then 
implies a bound on the value of the fundamental CP-odd parameter in the QCD action, 
6 < which is deemed to be unnaturally small. Since there is no good reason for this 

number to be so different from unity (he., a heretofore unknown symmetry in Nature), 
its minuteness requires “hne-tuning” of the action. This is often termed the Strong CP 
problem. A way around the hne-tuning is the so-called Peccei-Quinn mechanism based on 


a new (undiscovered) symmetry of Nature which requires that 6 be zero |11L 1121113L [1 
In this paper we present a calculation of in units of 6 within the framework of QCD 
with two havors of light quarks using the lattice regularization. A preliminary report on this 


work appears in the proceedings of the Lattice 2005 meeting m held at Trinity College, 
Dublin, and we note that while hnishing this work, a similar study, but in the quenched 
approximation, has appeared in M- 

As explained in Section II, the electric dipole moment is sensitive to the topology of the 
gauge held, or more specihcally, huctuations of topological charge; thus we focus mainly on 
calculations with dynamical, or sea, quarks. The two havor ensemble of lattice gauge helds 
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that we use was generated by the RIKEN BNL Columbia (RBC) collaboration. Details of 
these simulations are described in 17|. We hnd that a precise and accurate calculation re¬ 
quires ensembles with signihcantly longer evolutions {i.e., more independent conhgurations) 
than presently available; the topological charge has very long auto-correlations. We expect 
that longer evolutions will be available in the near future^. This situation is to be compared 
to the quenched, or zero flavor, case where topological charge can be evolved more efficiently. 
The topological charge susceptibility, however, does not vanish as the valence quark mass 
approaches the chiral limit, and as we show, neither does the electric dipole moment. This 
quenched pathology means can only be accurately calculated when the sea quarks are 
included 1^. Not surprisingly, this was found to be true in a recent work using the instanton 
liquid model |l^ where it was argued that the quenched chiral limit of is singular. The 
partially quenched limit rrivai —*■ 0, nisea fixed is also singular Q- 

Since topology is crucial in the calculation of djq, it may also be important to use lattice 
fermions that do not spoil certain topological relations to the gauge field with large lattice 
spacing errors. The axial anomaly in QCD relates the topological charge to the pseudoscalar 
density; a chiral rotation on the quark fields in the QCD action shifts the CP-odd 9 term 
between gluon and quark sectors. In order to realize this proper behavior, we use domain 
wall fermions which are chirally symmetric even when the lattice spacing is non-vanishing. 
Thus, this important continuum property of QCD is realized at non-zero lattice spacing, a 
feature that is absent for Wilson- and staggered- type fermions. 

This paper is organized as follows. Section |Tll describes the method to calculate ^at. 
Section IIIII gives details of the simulations, in Section US we present our results, and in 
Section El we summarize the present study and outline future calculations. 


^ The RBC and UKQCD collaborations are jointly beginning extensive simulations with 2-1-1 flavors of 
domain wall fermions. 
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II. GENERAL ANALYTIC FRAMEWORK 


A. Theoretical background 


We begin by considering the addition of a T- and P-odd (therefore CP-odd) term to the 
QCD Lagrangian (our conventions are detailed in Appendix A): 

SQCD,e = -9 j dt J d^x , ( 1 ) 

where G^y is the gluon held strength and the trace is over (suppressed) color indices. Such 
a term is allowed by the gauge, Lorentz, and discrete symmetries of QCD. It is easy to see 
that this so-called 9 term is odd under P and T transformations since 


e^^upaGP^G^^ ^ E-B. (2) 

0 is a fundamental, but unknown parameter of QCD. Remarkably, even though the r.h.s. 
of E^uafioomcao be woiften ae a fota, diveo.eace, . doee aof vaa^bQ and therefore has 
physical consequences, most notably CP violation in QCD. We return to this shortly. 

On the other hand, the QCD Lagrangian for massless fermions, 

^QCDJ = (3) 

is invariant under chiral transformations of the quark helds. 


V’ ^ (1+i075/2)V', 
i! - 0(1 + ^ 075 / 2 ), 


but the measure of the path intergral is not 


m, 


V'lpV'ij) —> VipV'ip exp 


d^x^^ii[e^^p^Gp,Gp„] 


( 4 ) 

( 5 ) 

( 6 ) 


which gives rise to the Adler-Bell-Jackiw anomaly. It is well known that this axial anomaly 
induces observable effects even in the CP even case, the mass of the r]' and the (relatively) 
large decay rate for 7 r° —> 77 to name just two. In this work we are interested in CP-violating 
effects, ones that vanish when the 9 term is absent from the Lagrangian. 

Choosing 0 = — 6 *, the 9 term can be rotated away, or canceled in the action. Recently, 
Creutz has proposed a scenario in the one-flavor theory where the 9 term can not be removed. 
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even in the massless limit j2j, |23| . Since we deal with at least two flavors of quarks, we will 
not consider this possibility further. 

If all the quarks are massive, the chiral rotation generates another term in the action that 
can not be canceled by further held re-definitions. 




(7) 


which is also P- and T- odd. Thus the CP-violating term in the QCD Lagrangian can be 
moved between the gauge and fermion sectors, but it can not be eliminated. 

Even though it cannot be eliminated, the 9 term can be written as a total derivative, as 
mentioned above. Still, as is well known, it does not vanish. For QCD 0 

/ = <3, ( 8 ) 

where Q is the integral topological charge {Q = 0, ±1, ± 2 ,...) . Thus the 6 term produces 
physical effects, like an electric dipole moment for the neutron. 

Theoretical calculations naturally yield in units of the unknown fundamental parame¬ 
ter 6. Thus to translate the current experimental bound to a constraint on 6, or to determine 
6 should a non-zero value of be found through experiment, requires evaluation of nucleon 
matrix elements. The lattice regularization of QCD provides a first-principles technique for 
such calculations which we describe after discussing the chiral limit of d^. 


B. Taking the chiral limit 

Consider the QCD partition function for Nf flavors of massive, degenerate, quarks after 
integrating over the Grassman quark fields, 

Z = j VAfj, det[.p(m, A^) +i9m'^^]^^ (9) 

where Afj, is the gluon held, lp{m, Afj) is a general covariant Dirac operator for a single 
quark flavor with mass m as may be found in the continuum or on the lattice. The choice 
of degenerate quarks is made for convenience of notation with no loss of generality in the 
following. Factoring out det lp{m, *^ 34 )) fhe CP-even part of fermion action, and assuming 6 
is small, 

det [Ip{m, A/j) + iOm'j^] = det[.^m, A^)] [1 -|- i9m tr( 75 .^m)“^) ] -|- 0{9‘^). (10) 
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Next, using the spectral decomposition of the inverse Dirac operator, 


JjXm, A^) ^ = ^- 


|A)(A| ^ 

iA + m’ 


( 11 ) 


where A is an eigenvalue and |A) an eigenvector of Ip, and the Atiyah-Singer index theorem, 
we hnd 


tr [75^‘(m/)] = ^ 

Q 


( 12 ) 

(13) 


m 

where n± are the number of right- and left-handed chiral zero-modes of p for a given gauge 
held conhguration with topological charge Q. In the limit m —>■ 0, the 9 term does not vanish 
from the action since the factors of m cancel, contradicting our above argument derived in 
the explicit m = 0 limit in Section III Al Similarly, it is not obvious that the held strength 9 
term, i9Q, vanishes as m ^ 0. 

The seeming contradiction is easily resolved by examining the role of the usual (CP-even) 
fermionic action. 


det lp{m, AP)^^ = 11 ^ {i\j + . 


(14) 


As m —>• 0 , Q 7 ^ 0 conhgurations are suppressed since they support exact zero-modes of p 
with zero eigenvalue. In other words, in the chiral limit the Q 7 ^ 0 conhgurations represent 
a set of measure zero, and the distribution of topological charge becomes a delta-function, 
5{Q)i with zero width, (Q^) = 0, so the 9 term ehectively vanishes. 

The quenched approximation of Equation UHl detppn, A^) = 1, still allows CP-violating 
physics since the pseudoscalar density term in the action (or equivalent y, the CP odd held- 
strength term) is not discarded (the same conclusion was reached in 2 ^ through a diherent 
line of reasoning). However, in light of the arguments just made, the mass dependence of 
any observable depending on 9 will be completely wrong, and one should expect signihcant 
systematic errors as a result. Indeed, the topological charge susceptibility, (Q^)/IA, which 
is closely related to as we have just seen, is well known to be non-vanishing in the pure 
gauge theory j^. One must also be careful in partially quenched theories {rrisea 7 ^ ^vai)'-, 


in a recent paper it was shown that the leading valence quark mass contribution to dj^ 
in partially quenched chiral perturbation theory is proportional to m^ea log Thus the 
limit rriyai —> 0 , rrisea fixed, is singular. 
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C. Computational Methodology 


The calculation of the dipole moments centers on the form factors that parameterize the 
matrix element of the electromagnetic current between nucleon states in the 6 vacuum, 


{p',s'\J^^\p,s)e = Us'{p')T^{q^)us{]^ 

+i 


F2{q^] 


(15) 

(16) 


2m 


Fsiq^) 


+(T^''qvl 


where 


= -u'y^u 
o 


2m 


-d'y^d. 


( 17 ) 


q = p' — p is the space-like momentum (g^ < 0 ) transfered by the external photon, s (s') 
the spin of the incoming (outgoing) nucleon, m is the nucleon mass, = z /2 [ 7 ^, 7 ^], and 
Us{p), Us{p) are Dirac spinors. 

The four terms on the right-hand side of Equation are the most general set consistent 
with the Lorentz, gauge, and CPT symmetries of QCD. The insertion of probes the 
electromagnetic structure of the nucleon; for g^ ^ 0 it is easy to show that Ei( 0 ) is the 
electric charge of the nucleon in units of e (-fl for the proton, 0 for the neutron), ^ 2 ( 6 ) is 
the anomalous part of the magnetic moment. Fa is the anapole moment, and ^ 3 ( 6 ) gives 
the electric dipole moment. The last two vanish when 6 —>■ 0. 

Later it will be useful to separate into its iso-scalar and iso-vector components. 




( 18 ) 


= u'y^u — d'f^d 

( 19 ) 


= u'y^u -|- d'y^d. 

( 20 ) 


D. Calculating dipole moments on the lattice 

Instead of directly computing matrix elements, lattice calculations proceed by studying 
correlation functions in Euclidean space-time. The desired S-matrix element is obtained 
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through the usual LSZ reduction formula, but in Euclidean space-time, usually by relying 
on the exponential dominance of the ground state, though in principle excited state con¬ 
tributions can also be obtained. For the case at hand, we study a three-point correlation 
function where a nucleon with spatial momentum p is created at time 0 by the interpolating 
held xjv(0,p), the current is inserted at time t, and then the nucleon state is annihilated at 
time t'. 

G>^{t',t) = J^{t,q)xNi^,P))- (21) 

Inserting a complete set of relativistically normalized states between each interpolating held 
and the current and translating all helds to eqnal times, we obtain 

= G^{q)x + (23) 

where .. ” denote excited state contributions which we ignore. Note, the correlation fnnc- 
tion contains the desired S-matrix element, with no need to analytically continne back to 
Minkowski space-time. For convenience, we separate the correlation fnnction into two parts, 
G^{q) which is a Dirac matrix, and f{t,f,E,E') which collects all the kinematical factors, 
normalization of states, and time dependence of the correlation function. Color indices have 
been suppressed. The interpolating held xn is the conventional one used in most lattice 
simulations, e.g., for the proton 

Xp = (^abc [u^Cx^db] Uc, (24) 

with a, b, and c color indices. 

The states are normalized conventionally, 

(0|xjvb,s) = \^Us{p), (25) 

and using the spinor relation^ 

'^Us{p)Us{p) = E{p)'y^ - ix ■ p + m (26) 


^ In this section we use Euclidean space conventions for the gamma matrices. See Appendix A for details. 
E{p) appears without a factor of i since we work in Euclidean time and momentum 3-space, that is, our 
external states are on-shell as they must be. 



and the projector 


(27) 

setting p' = 0 and 6* = 0, we find the magnetic form factor 

tiV^yG^{q^) = p,m(Fi(g2)+F2(g2)) (28) 

irV^VG^iq^) = -p.m(Fi(g2) + F2(g2)) (29) 

GM{q^) = F,{q^) + F^{q^). (30) 


Similarly, 


trP* G\q^) 


pi 

GEi.f) 


m{E + m) {F,(q^) + , 

11+7' 


4 2 ’ 

[Imp 


( 31 ) 

( 32 ) 

( 33 ) 


for the electric form factor GE{q‘^)- Thronghont this paper we inclnde the factor (1 + 7*)/2 
in projectors to yield the positive parity state (nentron or proton in the CP-even vacnnm) 
(see, e.g. Q)- 

To determine the desired moment, or form factor, the factor f{t,t',E,E') appearing in 
Eqnation 1221 must be removed from the correlation function. This is most easily done by 
taking a ratio with another suitably chosen three-point function. For example, taking the 
ratio of Equation |2H1 with Equation [^yields the magnetic dipole moment of the nucleon in 
the limit 0. 


lim — 

Py 


trP^^G|>jv(t,P,E,p) 

tiV^ Gp{t,t',E,p) 


lim - 

<j^o E + m 


Fijq") + 0 (<?") 

C^Piq^) 


1 tiV^yG^^j^jq^) 

Py tr'P*Gp(g2) ■ 

1 F,[q^) + F,[q^) ^ 

E + m G^p{q^) 


1 

2m 


1 + 0‘^^,P 


( 34 ) 

( 35 ) 

( 36 ) 


where we have used Ei(0) = 1 for the proton and 0 for the neutron, and = + 2 ( 0 ) is 
the anomalous part of the moment. P and N denote proton and neutron, respectively, and 
the denominator is always evaluated for the proton. Because we take ratios corresponding 
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to different components of the electromagnetic current, the hnite renormalization constant 
associated with the local lattice current^ drops out and need not be calculated. 


E. CP violating vacuum, 0/0 


In this section we consider the case 0/0. First, we must explain a somewhat subtle issue 
concerning mixing of the magnetic and electric dipole moment terms in correlation functions 


arising from the physical mixing of the 0 = 0 eigens 
is, the neutron mixes with P-odd states when 0/0 


;ates in the CP broken vacuum. That 
^ |27|^. This physical mixing of states 


gives rise to an rtnphysical mixing of the electric and magnetic dipole moment form factors in 
correlation functions like those given in Equation |221 Generally, the mixing can be written 
as a Dirac spinor with phase since 75 ^ 5 (p) = ns(p), he., 75 takes a spinor of a given 
parity into the other. So, instead of the spinor relation one obtains 


= E(p)7t - iy ■ p + 

s,s' 

^ E{p)'yt -i^ ■ p + m(l + 2*075), 
Us,eil^ = 


(37) 


(38) 


In the second line we have assumed that a -C 1. Using Ug^ip) in (1^ instead of Ugip) and 
proceeding as in the previous section, we obtain 

trp>^3/(/^(g2) ^ pym{Fi{q^) + F2{q^)) +VxVz{^F^ - 2irri^FA) + ^PxPzF2 + 0(0^) (39) 
trV^VGyiq^) = -p,m(Fi(g2) + _ 2tm‘^FA) + ^PyPzF2 + (P(02)(4O) 


Pz 


Pz 


tiV^y G/g/ = am{F - m)Fi + a{m{F - m) + ^)F 2 + + 0(0/ 


(41) 


trP-^G/g/ = (^amF,iq^) + a^^^F2iq‘^) + ^^Fsiq^)^+Oie^). (42) 

The hrst two equations are given for completeness; they are not used in our analysis. The last 
two can be used to extract the electric dipole form factor F^{q^). In particular, taking the 


^ On the lattice only the point-split fnr rn nf the current is conserved. Here we use a local definition, 

^ In a preliminary report on this workjl we did not account for this mixing. We are grateful to M. 
Pospelov and S. Aoki for pointing this out to us. The mixing occurs because we work with correlation 
functions, not directly with matrix elements constructed from CP eigenstates. 
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ratio of Equation 021 with the proton electric form factor correlation function, Equation 
we arrive at 


1 ^ 1 tiV-y 

ipz tTV^G^p{t,t',E,p} ipz tr'P*G'p(g2) 

_ amFi(g^) + aM^F2{q^) + 
m{E + m) G^p{q^) 


(43) 

(44) 


Subtracting the Ei and E 2 terms and taking the limit —>■ 0 yields the electric dipole 

moment. 


Fsiq^) 


1 tTV^yG%{t,f,E,^ amEi{q^) + a^^E2iq^^ 


2mGP{q^) [iPz tYV^G^p{t,f, E,p) m{E + m)GP{q‘^) 


= lim 


1 trP“’»G^((,(',£),3 amFj(q^) + aSi^F2(q 


( 45 ) 

( 46 ) 


2m F-o^ip^ tiP‘G‘p(t,t',E,p) m(E + m)GP 

The value of the mixing angle a is most easily calculated from the ratio of the zero 


momentum two-point functions 


27|. 




= Z 


NO 




2m^0 
Ws,e(0)h5',e (0) 
2m^e 


( 47 ) 






where, as usual, “... ” denotes excited state contributions. Using the spinor relation (EH) 
and appropriate projectors. 


(Xve(^)x]v9(0))e - Znc (48) 

(Xve(^)x]v9(0))e ~ iZNae~^^\ (49) 


to lowest order in 9. Note that Zj^e = Zn + O^O"^) and = rriN + OiO"^), and as 9 is 
very small in Nature, we work only to lowest order. Of course, the right-hand side of the 
first equation is nothing but the usual ground state contribution to the nucleon two-point 
function computed in the CP even vacuum. 

Some final remarks are in order. CP symmetry for fermions is conventionally defined 
assuming a real fermion mass. This is the condition that gives the form factors in Equation 117)1 
their usual interpretations, in particular, that the electric dipole moment is related to F 3 . 
Thus, in order to use Equation one must work in the standard basis. On the other 
hand, following j^, one may include in the action the mass term, arising from 
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a chiral rotation on the quark helds through a particular choice of basis, in addition to 
the —iQQ term used here. There it is shown that depends only on the combination of 
(renormalized) parameters Q = 6 + 6', and mixing effects like those described above will 
differ in just the right way to ensure this is so. In other words, the chiral rotation affects 
the quark helds in the correlation function as well as in the action. It is only the relative 
strength, or difference (note the opposite signs of the two terms), of the two contributions 
that leads to physical effects. That the physical value of the CP violating parameter must 
be 6 is clear since through renormalization QCD with bare 6^0 will generate a 6' term, 
even if the bare value of 6' is set to zero in the action (and visa-versa). Thus, transforming 
back to the standard basis, one arrives at the combination of renormalized parameters 6 + 6' 
which couples to GG. In the rest of the paper we will not use the notation 6, though it 
should be understood that when referring to the parameter 6 this is what is meant. 


F. Computing with 0/0 


The 0/0 action, being complex, is difficult to simulate with conventional lattice meth¬ 
ods. However, this problem can be avoided by working in the small 6 limit, 

1 


{ 0)0 = 


Z(6) 

1 


m 

= (O) - i6{QO) 


- i6Q)0 


(50) 

(51) 


where (P is a generic operator functional of the fields. Note { 0)0 becomes an expectation 
value in the CP-even vacuum, the CP-odd part weighted over topological sectors® 

{QO) = Y,P{QGQAO)u, (52) 

V 

where P{Q) is the probability that the gauge field configuration has charge Q. As before, 
the electric dipole moment, or any CP-odd observable, is seen to be closely related to the 
topological charge, and we expect that any such observable should vanish as {Q'^)/V —> 0. 
In this was shown explicitly for the large N limit. 


® In [23 it has been proposed to use the pseudo-scalar density as a weight instead. For chirally symmetric 
lattice fermions that have an index, this is equivalent to weighting with Q. If chiral symmetry is broken, 
then the two methods should agree in the limit a —> 0. 
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Chiral perturbation theory shows that ~ ml logm^ j6| and {Q’^)/V ~ ml [^, so each 
vanishes in the chiral limit, as expected. We will need the formula for the susceptibility, 

m Prnl 


1 / 


(53) 


later to compare to the lattice results, f ^ = 130.7 MeV is the pion decay constant to 

lowest order in chiral perturbation theory. 

Finally, the mixing angle a must also vanish as ml —> 0 since it is proportional to 6] 
this will happen as {Q'^)/V —>• 0. It bears repeating that in the quenched case {Q"^)/V is 
independent of the quark mass, implying that djsr and a do not vanish in the chiral limit. 

III. LATTICE DETAILS 


In this pilot study, the dipole moments of the neutron and proton are computed mainly 
using an ensemble oi Nf = 2 flavor QCD gluon conhgurations generated by the RBC collab¬ 
oration usin g do main wall fermions (DWF) and the doubly blocked Wilson (DBW2) gauge 
action. See [l2| for details and other basic physics results. Conhgurations were generated 
for three values of the sea quark mass, msea = 0.02, 0.03, and 0.04. In physical units, this 
range corresponds roughly to mg/2 < mgea ^ mg, where mg is the strange quark mass at scale 
H ^ 2 GeV. This study has focused on the = 0.03 and 0.04 ensembles, and we obtain 
results only for equal valence and sea quark masses, mf = m^ai = mg^a- 

Two- and three-point nucleon correlation functions have been computed on 440(460) 
lattices with sea quark mass mf = 0.03 (0.04), lattice volume 16^ x 32 sites, bare gauge 
coupling (3 = 0.80 (inverse lattice spacing a~^ ~ 1.7 GeV), hfth dimension size Lg = 12 , 
and domain wall height M 5 = 1.8. The lattices were generated with the exact 0 algorithm. 
Observables were computed on every 20 th trajectory except for the hrst 1650 trajectories 
at mgea = 0.04 where the separation was 25 lattices. 

Besides the dynamical calculations, we have computed two-point functions on 297 
quenched DBW2 [3 = 0.8^a“^ ~ 1.3 GeV) conhgurations, also generated and studied 
by the RBG collaboration j^. Here, m^ai = 0.05 (roughly mg). 

Throughout, periodic in space and anti-periodic in time boundary conditions are applied 
to the fermion helds. The gauge helds were generated with periodic boundary conditions. 


For the three-point functions, we use the sequential source method described in 


321, and 
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the source and sink are both Gaussian smeared, 



( 54 ) 


N is the number of times the smearing kernel acts on x, and a; is the width of the Gaussian 
that results in the limit N —> cx). We took iV = 30 and ut = 4.35, optimal parameters 


Q 


for quenched Wilson fermions and Wilson gluons at a ^ ^ 2 GeV 


. As it turned out, 


these were not optimal in the Nf = 2 DWF case, though they yield satisfactory results 
when compared to simple wall or point source interpolating helds. On the other hand, 
these parameters worked exceptionally well for quenched DBW2 (3 = 0.87. The three-point 
correlation functions were computed for two source times, alternating on successive lattices, 
t = 0 and 15 with sink times t = 10 and 25, respectively. The electromagnetic current 
was then inserted in the correlation function between these hxed source and sink times. 
Gorrelation functions on successive lattices were blocked together to reduce auto-correlations, 
leading to 220(240) pseudo-independent measurements in each case for rrisea = 0.03 (0.04). 
This is roughly the same blocking factor used in jl^ . 

To save computer time, the three-point functions were calculated from the non-relativistic 
components of the Dirac spinor which is equivalent to using the projector (1 -|- 7*)/2 on the 
source. To calculate the GP odd piece of the two-point function, the full four component 
Dirac spinor is required since using only the non-relativistic components gives identically 
zero for the projected two-point function. 

Since the magnetic and electric dipole moment terms in Equation are proportional to 
q'^, to compute the dipole moments, the correlation functions must be calculated for g 7 ^ 0 
and then extrapolated to g = 0 . Attempts to calculate wi th = 0 , for example by taking 
a derivative with respect to g, do not work on the latticeFor simplicity, we take the 
out-going nucleon to be at rest, and the incoming nucleon to have spatial momenta p = —g. 
Because high momentum states are more noisy, we restrict ourselves to total momentum 
IpP with Pi = 27r(?7,j)/Lj and — 4- Momentum conservation is enforced by summing 

over the location of the center of the smeared sink and Fourier transforming with respect to 
the current insertion point. 
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IV. RESULTS 


We begin this section by investigating the topological charge on the ensemble of iVy = 2 
gauge conhgurations. Figure Q shows the simulation time history of Q] evidently there 
are long autocorrelations, a fact already noted in ^|. The lower panel corresponds to 
a quenched simulation (a“^ ~ 1.3 GeV) where Q fluctuates rapidly. Note the abscissa 
is different in the quenched case. The difference in fluctuations reflects the fact that the 
quenched lattices are separated by 1000 sweeps, whereas the dynamical ones are separated by 
only hve trajectories, owing to the signihcantly higher cost of the latter. In the former case, 
one sweep consists of one heat-bath plus four over-relaxed hits on each link of the lattice. 
In the latter, one trajectory = 50 steps of hybrid molecular dynamics evolution of each link 
plus one global Metropolis accept/reject step. We also emphasize that the suppression of 
tunneling between topological sectors is an algorithmic, not physics, problem which is much 
worse in the dynamical case due to the smooth hamiltonian evolution of the Monte-Carlo 
algorithm (see also j^. 1^ for earlier studies of this problem using staggered fermions). 
The method used to calculate Q uses APE smearing with coefficient 0.45 for twenty steps 
and an improved dehnition of the lattice held strength (see for details). An even better 
approach may be to use the overlap dehnition of the topological charge js?], |38[ , though the 
precise dehnition of Q is probably not the limiting factor. 

In Figure El the topological susceptiblity x is shown for both quenched and Nf = 2 cases, 
the former being plotted as a horizontal line since it does not depend on any sea quark mass 


(the Nf = 2 results were determined from the data in 17|, the quenched from 


H). 


X and 


;o the appropriate power to make each 
17\ {Nf = 2) and 311 (quenched). The 


tUtt are plotted in units of the Sommer scale, tq, 
dimensionless. The values for ro were taken from 
interesting feature to note is the signihcant decrease of the Nf = 2 value relative to the 
quenched one. While there may be some sea quark mass dependence, x levels oh between 
'msea = 0.03 and 0.04. In addition, the statistical errors shown in the hgure were estimated 
by blocking the data in groups of 50 trajectories (10 lattices) and treating the blocks as 
independent while Figure ^ indicates the topological charge has autocorrelations on longer 
scales. Also shown in Figure El is the prediction from lowest order chiral perturbation theory 
(Equation ESI), is comforting that this lowest order prediction is consistent with the 
Nf = 2 lattice calculation, but because of the caveats just mentioned, the agreement is 
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not yet significant. Given the close relation between y and the quark mass dependence of 
the electric dipole moment, it does not appear promising that the mass dependence of 
can be accurately determined from these ensembles; (much) longer evolutions are required. 
Nevertheless, the frif = 0.03 and 0.04 calculations may give a relatively good estimate of the 
magnitude of in QCD where the lightest quark mass is about m^. From Figure |2l this is 
almost surely not true for the quenched case. 

Next we discuss the CP even and odd parts of the two-point function (Equa¬ 
tions EH and EHl)- Again, working to lowest order in 9 by weighting expectation values 
with iQ in the latter case, the masses and Z factors obtained from each must be equal. To 
reduce statistical errors, we average the forward and backward in time parts of the nucleon 
propagator and over source time slices 0 and 15 for the mgea = 0.03 calculation®. For the 
usual 6 = 0 propagator, the former means averaging positive and negative parity states 
(particle and anti-particle). For 6^0, the particle and anti-particle states have the same 
CP-odd part containing both parities (c./.. Equation l37jl . Thus, we £t to 

Geven(t,l^ = (55) 

G,d^{t,p) = A (e-^® * - , (56) 

for the former and latter, respectively, in the range 7 < t < 12 to avoid excited state 
contamination. Ignoring the excited state contributions is justihed by the acceptable y^/dof 
of the single particle hts (Tables HI and lflf|l . For the CP odd case, the y^/dof is a bit large in 
some cases (Tables HU and I1V|1 , but likely for different reasons that are explained below. The 
average over forward and backward propagating states is equivalent to performing a time- 
reversal transformation on the correlation function which, in turn, is equivalent to averaging 
over time-reversed gluon configurations. This last step flips the sign of the topological charge 
on the underlying gluon configuration (recall that the 6 term is odd under time-reversal). 
Thus, performing the average of forward and backward correlation functions has the same 
effect as exactly symmetrizing the topological charge distribution of the ensemble. 

TableHlandHTlsummarize the hts to the two-point function for the Nf = 2 case. For rUsea = 
0.04 the measured values of E{p) are compared to the continuum relativistic dispersion 

® We omitted the negative parity states from our earlier discussion and Equation El for clarity. The 
backward propagating, negative parity, anti-particle state appears because of the anti-periodic boundary 
condition in time (see 0 ). 
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relation, 


E{p) = + (57) 

in Figure ini with pi = 27m/Li and sin(j9j) (n = 0, ±1, ±2,±Lj — 1), the latter being 
the exact lattice momentum for a free lattice fermion with nearest neighbor action. The 
agreement is satisfactory for small |p|, indicating lattice artifacts are small in this case. 
E{p) and mAr in (jKTjl are taken from the CP even part of the correlation function. Thus, 
in the following we simply use pi = 2Tm/Li for the momentum. Since this would lead to 
large O(a^) errors for large \p\, we restrict our analysis to the four non-zero lowest values 
admitted on our lattice, n = (1,0,0), (1,1,0), (1,1,1), and (2,0,0), and permutations. Since 
the larger momentum correlation functions are considerably more noisy anyway and would 
suffer large 0{a‘^) errors with either choice, this is not a cause for concern. Finally, we 
note that the values of the nucleon masses give n in Table HI differ by about two (statistical) 
standard deviations from those reported in |l2|. Different fit ranges and sources were used 
(8 < t < 16, wall source in jl^), and our statistical errors are about two to three times 
smaller, reflected by our increased statistics. 

In Figure E] we show the htted nucleon mass, mgea = 0.03, versus the minimum time slice 
used in the £t. Values of mAr for both CP even and odd parts of the two-point function are 
shown. For tmin > 5, the masses are constant within statistical errors and agree with each 
other, as they should to this order in 6 for the CP odd part of the correlation function. Of 
course, the errors on the masses from the CP odd part are significantly larger, due to the 
topological charge re-weighting procedure. A naive ratio of the two-point functions which 
gives the mixing angle a is also shown in Figure 0] (lower panel). Again, for f > 5 a suitable 
constant plateau is evident from which we extract the average value a(0.03) = 0.16(2) 
(6 < t < 9). 

In Figure El we show the results of a similar analysis, but for m^ea = 0.04. This time, in 
the whole range of tmin the masses clearly disagree outside statistical errors; the difference 
in central values is roughly ten percent. The value of y^/dof for the CP-odd case is roughly 
two, while in the CP-even case it is less than one (Tables H] and EJ. A naive ratio of the 
two-point functions, which would give the mixing angle a if the masses were equal, is also 
shown in Figure El (lower panel). Although a plateau appears at small tmin, the ratio appears 
to decrease approximately linearly with t in the region where the masses are constant but 
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unequal, as expected. As mentioned above, for the nisea = 0.03 case, we were able to average 
over two source times, whereas for nisea = 0.04, we only calculated the full four component 
two point function needed for the CP odd part from one source time. This suggests that the 
different behavior between the two cases may be due to an improved overlap of the nucleon 
helds with the local charge in the former case and that even more source times would improve 
the CP odd signal. 

To further study these effects, we have calculated the same two-point functions on a 
quenched ensemble of lattices. Unfortunately, it was only after the nigea = 0.04 and quenched 
calculations that we realized the importance of having more than one source for the CP odd 
part of the correlator. As mentioned already, the topological charge distribution on this 
ensemble is expected to be correct in quenched QCD because many more Monte-Carlo 
updates have been performed between measurements. The masses obtained from hts like 
those in the dynamical case are shown in Figure Inland given in Tables EH and IIVI Note that 
the plateaus are quite good in this case, suggesting the Gaussian smearing parameters are 
near optimal. The masses agree within statistical errors for t > 5 and the difference of the 
central values is less than hve percent, so now the naive ratio provides a relatively accurate 
value of the mixing angle, a = 0.214(32), where the error is statistical only, and we have 
averaged over the range 5 < f < 10. All of the hts for both CP-even and odd parts of the 
correlation function have y^/dof < 1. While the quenched result is clearly an improvement 
over the two havor one. Figure IHl suggests that an even more accurate sampling of Q is 
desirable. In Figure E E{p) computed from the CP even part of the correlation function is 
shown. 

The agreement of the masses between the CP even and odd parts of the correlation 
function is a simple, but non-trivial, check that the distribution of Q is correct. Summarizing 
our initial studies, the CP odd signal may be improved signihcantly by using more sources 
for the correlator and by increasing the sampling of global topological charge. The former 
improves the overlap of the nucleon with the local charge fluctuations. 

One way to try to salvage the nigea = 0.04 calculation is to £t each correlator separately, 
extract the coefficient of the ground state exponential from the large time region and then 
take the ratio of these coefficients to determine a. One can take the mass in the CP odd 
case as a free parameter or £x it to the correct value from the CP even case. Though 
this is more correct than just taking the ratio of correlation functions and picking out the 
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incorrect plateau, it will still yield a value of a with some significant systematic error. After 
all, the htted masses differ by about ten percent. Likewise, we may anticipate the value of 
a to be incorrect by this amount. Being a bit more systematic, let us say the correlation 
function itself has been determined close its actual value. Then we make a ten percent error 
in the amplitude (and therefore a) since the htted mass is ten percent too low. Now, a 
reasonable guess may be that the correlation function is actually determined to roughly ten 
percent of its correct value. Taking both factors into account, we arrive at a systematic 
uncertainty in a of about 10-20%. Following this procedure gives a = 0.07(2) and 0.16(2) 
for rrisea = 0.04 and 0.03, respectively (statistical error only). Here, ttin in the CP odd part 
of the correlation function has been hxed to the CP even value, and a is from the ht with 
tmin = 7. The agreement with the simple ratio method in the case rrisea = 0.03 is a nice 
check of both methods. If mAr is left as a free parameter, the resulting value of a is about 
50% lower (Figure |H1) for rrisea = 0.04, but the same at 0.03 with larger error, a = 0.18(4). 
So while there is likely signihcant uncertainty in the value of a at rrisea = 0.04, the results 
for the lighter mass appear stable and satisfactory 

The ratios of the three-point correlation functions given in Equations El and |131 are 
shown in Figures IHl - El for each value of and both sea quark masses. For magnetic 
form factors we average results using Equations EH and EH Figure El shows the ratio 

using the value of a determined above. Despite the flat plateaus shown 
in the hgures, some excited state contamination may still be present at small and large times, 
given the htted masses in Figure El For some cases, involving larger momentum transfer, the 
plateaus show an oscillation, presumably due to insufficient statistics, or possibly excited 
state contamination. Averages over time slices 4 < t < 6 are also given in Table E The 
ratios Fi^ 2 {q^)/{q^) and F^{q^)/{2mGP{q^)) are summarized in Table IVTl The value 
of the Fi ratio for the neutron approaches zero, as required by electric charge conservation. 
For the proton, the ratio trivially approaches one since it goes to Fi(0)/Fi(0) (Equations EH 
and ESI), but at the least serves as a check on our evaluation of the three-point functions. 

^ While finishing this manuscript, we noticed a discrepancy between our value for a in the quenched case 
and that reported in 0 Our definition of the mixing angle differs from theirs by a factor of 1/2 (compare 
Equation E3 with Equation 20 in |lfilli. However, the quoted value in [7^ is numerically equal to ours, 
within statistical errors. Unfortunately, communication with the authors of [7^ and further checks of our 
definitions, code, and analysis of the two point functions has not yet resolved this discrepancy. 


19 



As —> 0, the F 2 {q^) ratios yield the anomalous magnetic moments of the nucleons. For 
each value of the magnitudes for the neutron and proton are equal within errors; this 
should be true for the iso-vector contributions, assuming iso-spin is not broken which is true 
in our calculation. Evidently the iso-scalar contribution from the connected diagrams is zero, 
or smaller than our statistical errors; we have not included the disconnected valence quark 
loop diagrams in the three-point functions which contribute only to the matrix element 
of the iso-scalar piece of the electromagnetic current. The values at the lowest value of 

are not far off from the well known experimentally measured values = 1-79 and 
= —1.91 1^. To show the momentum dependence of these ratios is mild and to 
compare to experiment j^, we have also plotted the ratio of the magnetic form factors to 
the electric form factor of the proton, GM{q^)/{q^) ^ as well as the dipole moments in 
Figure fT7l Note, the quark mass dependence of this ratio is also small, mainly showing up as 
shift in for larger values of q^. In fact, it is interesting that the last two mgea = 0.03 points 
seem to smoothly £11 in the large gap between the last two irtgea = 0.04 points. By comparing 
to experiment, one sees that lattice artifacts are becoming significant for g^> 1.2 (GeV^). 
The agreement with the experimental form factor ratio for the proton at smaller values 
of q^ is quite satisfactory, in magnitude and q^ dependence, but may be fortuitous since 
our calculation does not include electromagnetic effects or disconnected valence quark loop 
contributions, is done at relatively heavy quark mass, and we have not taken the continuum 
or infinite volume limits. 

Finally, we turn to the CP odd form factor F^{q^). The 0 limit of the F‘i{q^) ratio 
yields |dAf|- From the last two columns in Table lyTl it can be seen that our value for the 
form factor is consistent with zero within errors for both quark masses and both components 
of the electromagnet current, except for the second lowest value of g^ for each quark mass, 
where the central value is roughly two standard deviations from zero ( 7 * component only). 
The statistical errors are much larger for the z component of the electromagnetic current 
(Equation HU) compared to the t component (Equation |121), hence we discuss only the latter 
from now on. The statistical errors for the lighter quark mass are about twice the size for 
the heavier. For the three point correlation functions, we have averaged over two source-sink 
time slice combinations, with sources separated by 16 time slices, or one-half the lattice in 
the time direction. Thus the increased error at light quark mass may reflect the (usual) 
greater fluctuations in the quark propagator as the quark mass is reduced. We are also 
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mindful that lowering the quark mass suppresses the evolution of the topological charge and 
charge density and likely plays a role as well. Again we emphasize that future calculations 
may beneht (in terms of smaller statistical errors) from improving the correlation of the 
nucleon helds and the local charge density by using many source-sink combinations instead 
of the two used here. Indeed, the statistical error on the ratio already decreased by about a 
factor of two by using two sources. 

Given the large statistical uncertainties, we simply estimate a bound on from the error 
on the F^{q^) ratio evaluated at the lowest value of momentum transfer for mgea = 0.03. 
Though the error on the values at rusea = 0.04 are a factor of two smaller, we do not use 
these values because of the possibly large systematic error in the value of a extracted from 
the two point correlation function. We hnd 


d 


lat 

N 




^ FsjO) 

2m 

F3(0.399) 

^ 2mGij(0.399) 
< 0 . 2 . 


( 68 ) 

(59) 

( 60 ) 


The mild dependence for the Fi and F 2 ratios leads us to believe that this is not a terrible 
approximation. Our conventions which are the same as those in j^, 0| have lead to a negative 
central value of d^/O at msea = 0.03 and positive at 0.04. Of course, both are consistent 
with zero within errors. We simply mention this to call attention to our careful treatment 
of the many sources of convention-based signs in this calculation (see section II and the 
appendix) and so others may compare their calculations to this and future ones. 
d}^^ is given in inverse units of aO. In physical units, the above bound reads 


dN ^ 0.02 6 e fm. 


(61) 


This upper bound is roughly an order of m^nitude larger than (but not inconsistent with) 


the central value computed from sum rules [2| and a factor of four times as large as the pion 


loop contribution p. Also, the above bound is roughly the same size as the result found in 


a recent quenched calculation [l^; however, see 
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V. CONCLUSIONS 


The dipole moments of the proton and neutron have been calculated using lattice QCD. 
In particular we have focussed on the electric dipole moment of the neutron. Using two flavor 
QCD, we obtained a rough bound, djsf < 0.02e6*fm, from the statistical error on the central 
value which was zero within this error. This bound is somewhat larger than previous model 


Q. 


and 


calculations, about the same magnitude as found in a recent quenched calculation 
is only a crude estimate given additional significant systematic uncertainties associated with 
the topological charge distribution and quark mass and momentum dependence of F^{q^). 

The ratio of magnetic to proton electric form factors were found to be in good agreement 
with experiment fFigure IT7j) . given the single lattice spacing, volume, and relatively heavy 
quark masses used in our calculation. This ratio exhibits only a mild dependence on the 
momentum transfer and quark mass. 

Because djq arises from the CP-odd term in the action, J d'^xGG, it is sensitive to the 
topological charge distribution. The method outlined here to calculate CP odd observables 
requires re-weighting correlation functions with topological charge that would otherwise 
vanish. This suggests that reducing the statistical errors in such calculations can be achieved 
by improved sampling of the topological charge, and also by increasing the number of quark 
sources so the overlap between the nucleons and the local topological charge density is 
enhanced. The former was shown to work in the quenched approximation, while evidence 
for the latter came from the rusea = 0.03 calculation of the CP odd part of the two point 
correlation function. 

As discussed in Section |TT1 in the quenched case d^ does not have the correct (physical) 
quark mass dependence. In the two flavor case, the quark mass dependence is correct, 
and djsi vanishes in the chiral limit from the presence of the CP even part of the fermion 
determinant. Future lattice calculations will approach the chiral limit, which will sujppress 
even further. Chiral perturbation theory predicts a leading m^logm^ term p|, but 


non¬ 


leading ml terms may also be important p. 


A recent calculation in partially 


quenched chiral perturbation theory m may help with the needed extrapolations. It is also 
interesting to note the prediction in that paper for the leading valence quark mass dependent 
term, m^ea logThus the limit m^ai —^ 0, mgea fixed, is singular. 

Nf = 2 -|- 1 flavor domain wall fermion calculations just begun jointly by the RBC and 
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UKQCD collaborations will attempt to address the two most pressing deficiencies of the 
present calculation, poor statistics for the topological charge and the quark mass depen¬ 
dence of d]sf. We plan to implement the two main lessons learned from this study in the new 
one, namely longer evolutions of the gauge fields and more quark sources to overlap with the 
local charge density. Both should improve the statistical errors on CP odd correlation func¬ 
tions. Finally, we will investigate the use of twisted boundary conditions to essentially 
eliminate, or at least drastically reduce, the extrapolation of the form factors to q'^ = 0. 
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APPENDIX A: CONVENTIONS 

In Minkowski space it is conventional to dehne the chiral basis with metric (signature: 
1,-1,-1,-1) as 


7° = 


0 1 
1 0 , 


where = 


r = 


0 1 
1 0 



a2 = 


i 0 


a3 = 


1 0 
0 -1 


are the Pauli matrices. 


In reality, our code uses slightly diherent conventions. With the replacement je = 
and 7 I; = 7°, 


— ?7 = 


—t'f = 


-I'j = 


7 = 


-7i 

7I 

-71 

7e, 


(Al) 

(A2) 

(A3) 

(A4) 


he., our dehnitions of and 7 ^ have the opposite sign of the usual ones (still, = 
= — 7 ^)- Of course, the results of Eqs. EE EE IIE and EE do not depend on the 
choice of metric or signs for the gamma matrices as long as the projectors are also modihed 
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accordingly. We have checked this explicitly by comparing results from our conventions to 
ones using the above conventional chiral basis. In the above the subscript “E” stands for 
Euclidean and is dropped in the main text. 

In Minkowski space the gluon action takes the form 

,2 

^ flV 

±U /I “ 

dt J d^xC 


c = - »^G-G„ 


^ = 


dt I d^x ( + e-^G^^'^G 

4 lOTT"' 




Continuing to Euclidean space, 


(AS) 

(A6) 

(A7) 


t - 

->• —IT 

(AS) 

GOi _ 


(A9) 

Goi - 

-iG^^ 

(AlO) 

Gij _ 


(All) 

n 

Kjr^j 

-G*h 

(A12) 


The continuation of the held strength term can be worked out from 


^ -iA^ 

G^u = ^iiuagG^^j 


(A13) 

(A14) 

(A15) 


and similarly for the (contra)covariant (dual)held strength. A^ is the four-vector potential 
and 60123 = ei 234 = +1- Then 


1 „2 

expfA ^ exp <; - / dr / d^x ( 

4 lOTT"' 


(A16) 

(A17) 


= expj-S'ij}. 

Thus, the sign of the 6 term does not change upon analytic continuation to Euclidean space 
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TABLE I: Masses and energies from fits to the CP even part of the nucleon two-point function 

{0 = 0). Results are averaged over all possible permutations of the lattice momentum given in 

the first column, including both positive and negative directions. The fit range is 7 < t < 12. 

fnsea = 0.03 (upper) and rusea = 0.04 (lower). 

P 

E{p) (error) 

X^/dof 

dof 

(0,0,0) 

0.8646 (53) 

0.74 

4 

(1,0,0) 

0.9453 (65) 

1.6 

4 

(1,1,0) 

1.031 (10) 

1.3 

4 

(1,1,1) 

1.104 (18) 

2.1 

4 

(2,0,0) 

1.140 (29) 

1.3 

4 

(0,0,0) 

0.9264 (54) 

0.31 

4 

(1,0,0) 

1.0021 (61) 

0.39 

4 

(1,1,0) 

1.0685 (88) 

0.25 

4 

(1,1,1) 

1.124 (16) 

0.45 

4 

(2,0,0) 

1.217 (30) 

0.67 

4 


TABLE II: The same as Tabled except for the CP odd part of the two-point function {0 ^ 0). 

Note, the fit range for momentum (1,1,1) and mgea = 0.03 begins at t = 6. 

P 

E{p) (error) 

xVdof 

dof 

(0,0,0) 

0.881 (34) 

0.50 

4 

(1,0,0) 

0.974 (65) 

0.77 

4 

(1,1,0) 

1.10 (18) 

2.27 

4 

(1,1,1) 

1.17 (42) 

0.91 

5 

(0,0,0) 

0.814 (36) 

1.85 

4 

(1,0,0) 

0.873 (54) 

1.90 

4 

(1,1,0) 

0.950 (86) 

2.48 

4 

(1,1,1) 

1.18(16) 

0.93 

4 
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TABLE III: The same as Table QJ except for the quenched ensemble. 


p 

E{p) (error) 

xVdof (dof) 

(0,0,0) 

1.0217 (44) 

0.46 (4) 

(1,0,0) 

1.0889 (52) 

0.53 (4) 

(1,1,0) 

1.1538 (70) 

0.24 (4) 

(1,1,1) 

1.2152 (110) 

0.14 (4) 


TABLE IV: The same as Tabled except for the quenched ensemble and 0 0. 

P 

E{p) (error) 

xVdof (dof) 

(0,0,0) 

0.996 (30) 

0.97 (4) 

(1,0,0) 

1.054 (31) 

0.36 (4) 

(1,1,0) 

1.123 (47) 

0.71 (4) 


TABLE V: Ratios of three-point functions given in Equations I1I4I and I4dl In the limit 0, these 
yield the dipole moments, except in the case of the electric dipole moment (last column) where the 
mixing with the Fi and F 2 terms has not been subtracted, mgea = 0.03 (upper) and mgea = 0.04 
(lower). 


(GeV2) 

proton (magnetic) 

neutron (magnetic) 

neutron (electric) 

0.399 

1.524 (54) 

-0.957 (37) 

0.35 (19) 

0.824 

1.624 (65) 

-1.030 (43) 

-0.11 (22) 

1.183 

1.751 (129) 

-1.075 (84) 

0.43 (36) 

1.363 

1.408 (176) 

-0.952 (131) 

0.30 (53) 

0.401 

1.438 (30) 

-0.913 (23) 

-0.045 (87) 

0.753 

1.451 (34) 

-0.936 (26) 

0.066 (96) 

1.044 

1.594 (68) 

-1.032 (50) 

-0.035 (160) 

1.538 

1.193 (104) 

-0.730 (69) 

-0.393 (249) 
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TABLE VI: Form factors normalized by the electric form factor of the proton, GEiq^) (Equa¬ 
tion 01. In the limit 0, the values in the F 2 columns yield the anomalous magnetic moments 
and in the last two columns, the electric dipole moment of the neutron using currents and projectors 
defined in Equations ETI and W]\ rusea = 0.03 (upper) and nisea = 0.04 (lower). 


(GeV2) 


F2{q^)/GE{q^) 

(Fsiq^)/^ 

m)/GE{q^) 


proton 

neutron 

proton 

neutron 

neutron 

0.399 

1.0784 (63) 

-0.034(10) 

1.680 (98) 

-1.698 (68) 

-0.04 (20) 

0.49 (45) 

0.824 

1.183 (15) 

-0.1261 (211) 

1.90 (12) 

-1.827 (82) 

0.45 (23) 

1.56 (73) 

1.183 

1.297 (40) 

-0.219 (42) 

2.15 (23) 

-1.90 (15) 

-0.08 (36) 

-0.74 (1.73) 

1.363 

1.251 (57) 

-0.176 (71) 

1.57 (31) 

-1.73 (25) 

0.02 (53) 

-0.39(73) 

0.401 

1.0695 (42) 

-0.045 (7) 

1.703 (59) 

-1.715 (46) 

0.087 (95) 

0.12(27) 

0.753 

1.1349 (89) 

-0.081 (13) 

1.760 (66) 

-1.785 (52) 

0.20(10) 

-0.18(41) 

1.044 

1.2181 (223) 

-0.102 (27) 

2.050 (129) 

-2.013 (101) 

0.12 (16) 

-1.29 (75) 

1.538 

1.2107 (414) 

-0.156 (45) 

1.345 (195) 

-1.409 (136) 

-0.29 (25) 

0.55 (38) 


29 









Trajectory 


FIG. 1; Topological charge, Q. For the Nf = 2 simulations, Q for every fifth trajectory is shown, 
while for the quenched case Q has been measured on lattices separated by 1000 sweeps. The Nf = 2 
plots are reproduced from [17 1. 
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FIG. 2: Topological charge susceptibility for Nf = 2 (filled circles) and quenched (solid line and 
horizontal dashed lines) simulations shown in Figure ^ The dashed line is the chiral perturbation 
theory prediction, Equation [221 with r^f evaluated from Results are given in terms of the 
Sommer scale, tq, for convenience. 
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■lpl 


FIG. 3: The nucleon energy from a fit to Eauationl55l \p\ = + n| where re* = 0, ±1, ±2. 

The dashed line is plotted from the continuum formula, E{p) = and the solid one is 

the same except pi is replaced by sin(pj). All quantities are shown in lattice units, nisea = 0.04, 
a-i Ri 1.7 GeV. 
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FIG. 4: Upper panel: nucleon mass from single particle fits to the CP-even (cirlcles) and odd 
(squares) parts of the two-point correlation function; tmin is the value of the smallest time slice 
used in the fit. Lower panel: mixing angle a from the simple ratio of the same CP-odd and even 
parts of the two-point correlation function, msea = 0.03. 
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FIG. 5: Upper panel: nucleon mass from single particle fits to the CP-even (cirlcles) and odd 
(squares) parts of the two-point correlation function; tmin is the value of the smallest time slice 
used in the fit. Lower panel: mixing angle a from the simple ratio of the same CP-odd and even 
parts of the two-point correlation function, rusea = 0.04. 
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FIG. 6; Same as Figure El but for the quenched simulation described in the text. 
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FIG. 7: Same as Figure 01 except for the quenched ensemble described in the text. 
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FIG. 8: The mixing angle a computed from fits of the CP even and odd parts of the nucleon 
two-point function to Equations |1H1 and I4h[ Results are shown for mjv fixed in the latter case to 
the value from the 0 = 0 part (squares) and for mjv a free parameter, rusea = 0.04. 
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FIG. 9; Ratio of three-point functions given in Equation[^that yields the magnetic dipole moment 
of the neutron in the limit —> 0. Plots are shown for units of lattice momenta, q = (1,0,0), 
(1,1,0), and (1,1,1), and permutations, in the lower, middle, and upper panels respectively, rngea = 
0.03. 
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FIG. 10; Same as in Figure El but for the proton. 
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FIG. 11: Same as in Figure El but nisea = 0.04. 
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FIG. 12: Same as in Figure IHl but for the proton, and rusea = 0.04. 
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FIG. 13: Same as in Figure IHl but for the electric dipole moment of the neutron iFnuatiou 1431) . 
The mixing with the Fi and F 2 terms has not been subtracted, rusea = 0.03. 
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FIG. 14: Same as in Figure IHl but for the electric dipole moment of the neutron (Equation |321) • 
The mixing with the Fi and F 2 terms has not been subtracted, rusea = 0.04. 
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FIG. 15: The ratio {F2,{q^)/2m)/GE{(f‘) (Equation ESJ . increases from bottom to top panel. 
In the limit —> 0 this ratio yields the electric dipole moment of the neutron. Also shown is the 
subtraction term (squares) in Equation rngea = 0.03 
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FIG. 16: Same as in Figure ITKl but rusea = 0.04. 
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FIG. 17: Ratios of the proton (squares) and neutron (circles) magnetic form factors to the proton 
electric form factor (Equation 03 times E + m). The limit ^ 0 yields the magnetic moments. 
The absolute value is plotted for the neutron for comparison, rn™ = 0.03 (filled symbols) and 
0.04 (open symbols). The diamonds are experimental data points |4fll| where we have added quoted 
statistical and systematic errors in quadrature. The triangles at = 0 are the experimentally 


measured magnetic moments m (the absolute value for the neutron (lower triangle) is shown). 
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